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ABSTRACT 


Several  experiments  have  suggested  that  the  streamwise  vortices,  with  their 
accompanying  low  momentum  streaks,  in  a  turbulent  boundary  layer  have  a  char¬ 
acteristic  spanwise  wavelength  of  approximately  =  100.  A  mechanism  is 
proposed  which  selects  a  comparable  spanwise  wavelength  and  produces  counter¬ 
rotating  streamwise  vortices  in  a  turbulent  boundary  layer.  Examining  the 
equations  which  describe  the  small  deviation  of  the  velocity  field  from  its 
time-average,  it  is  found  that  the  Benney-Gustavsson  resonance  (Studies  in 
Applied  Mathematics  1981)  occurs  with  such  a  boundary  layer  velocity  pro¬ 
file.  It  is  shown  that,  as  an  integral  part  of  this  resonance,  there  is  a 
mean  secondary  flow  which  has  a  spanwise  wavelength  X^  =  90  and  whose 
velocities  exhibit  a  counter-rotating  streamwise  vortex  structure. 
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1.  INTRODUCTIOM 

During  the  last  decade,  much  of  the  research  in  turbulent  boundary  layers  has 
been  concerned  with  the  coherent  eddy  structures  that  have  been  observed  near 
the  rigid  boundary.  The  quasi-deterministic,  randomly-located  sequence  of 
coherent  structures,  collectively  called  the  bursting  process,  is  believed  to 
play  a  dominant  role  in  the  production  and  maintenance  of  turbulent  flow  (for 
a  review,  see  Cantwell  [1]). 

One  important  aspect  of  the  bursting  process  is  believed,  by  many  investiga¬ 
tors,  to  be  the  counter-rotating  streamwise  vortex  structure  with  its  accom¬ 
panying  low-speed  streaks  (for  example,  see  Blackwelder  [2]).  From  the  many 
measurements  of  these  streaks,  it  has  been  found  that  the  mean  streak  spacing 
is  approximately  100  v/p/t  ',  where  v  is  the  kinematic  viscosity,  p  is  the 
density  and  is  the  shear  stress  at  the  wall.  From  flow  visualizations  [3], 
these  low-speed  streaks  are  usually  observed  to  end  by  being  lifted  av/ay  from 
the  wall.  As  they  are  lifted,  the  streaks  start  to  oscillate  with  increasing 
amplitude  leading  to  what  is  called  breakdown.  Corino  and  Brodkey  [4]  showed 
that,  soon  after  this  breakdown,  a  large  scale  motion  emanating  from  the  outer 
flowfield  approaches  the  wall  and  "cleans"  the  entire  region  of  the  chaotic 
motion.  This  phase  of  the  bursting  process  has  been  called  the  sweep. 

In  experiments  using  hot-wire  anemometers,  the  bursting  process  is  most  easily 
characterized  by  a  sharp  acceleration  or  the  streamwise  velocity  as  reported 
by  Wallace  et  al.  [5]  and  Blackwelder  and  Kaplan  [6].  To  single  out  the 
bursting  process  from  the  chaotic  turbulent  motions,  Blackwelder  and  Kaplan 
used  a  burst  detection  technique  which  relies  on  this  feature  of  the  bursting 
process.  From  conditional  ensemble  averages  of  data  taken  in  such  a  way, 
Blackwelder  suggested  that  the  sharp  acceleration  is  the  transition  between 
the  low-speed  streak  and  the  sweep. 

Recently,  Blackwelder  emphasized  the  similarity  between  the  turbulent  boundary 
layer  bursting  process  and  the  transition  problem  [7].  In  particular,  he 
found  that  the  conditionally  averaged  streamwise  velocity  profile  is  highly 
inflectional  at  and  before  the  burst  detection.  He  further  speculated  that  a 
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localized  shear  layer  instability  might  play  a  role  as  in  the  transition  prob¬ 
lem.  The  spatial  scale  of  the  oscillation  of  the  lifted  streak  observed  in 
flow  visualizations  was  noted  to  be  consistent  with  such  an  interpretation. 
In  any  case,  in  Blackwel der ' s  model,  the  counter-rotating  streamwise  vortex 
pair  plays  several  important  roles.  The  low-speed  streaks  are  attributed  to 
the  vortex  pair '  s' ' pumpi ng '  action  of  the  low  speed  fluid  av;ay  from  the  v;all. 
In  turn,  the  low-speed  region  is  responsible  (at  least  partially)  for  the 
inflectional  velocity  profile  which  eventually  leads  to  breakdown  through  a 
shear-layer  instability  mechanism. 

The  aforementioned  similarity  between  the  bursting  process  and  the  transition 
problem  has  suggested  to  some  that  the  streamwise  vortices  observed  in  either 
case  might  be  due  to  the  same  mechanism.  In  this  vein.  Coles  [8]  speculated 
that  the  streamwise  vortices  in  a  turbulent  boundary  layer  might  be  the  result 
of  a  Taylor-Gortler  type  instability  in  which  the  concave  flow  is  due  to  the 
large  scale  motion  in  the  outer  flowfield.  This  idea  predicts  a  spanwise 
wavenumber  for  the  vortex  structure  consistent  with  experiment  [2].  However, 
it  appears  difficult  to  justify  the  extention  of  the  steady  state  Taylor- 
Gortler  stability  analysis  to  the  unsteady  flowfield  near  the  flat  wall. 

Benney  and  Lin  [9,10]  proposed  another  mechanism  which  attributes  the  vortex 
structure  in  the  transitional  flow  to  the  secondary  mean  flow  induced  by  the 
nonlinear  interaction  of  a  two  and  a  three-dimensional  wave.  Although  this 
mechanism  reproduced  the  essential  features  of  Klebanoff's  experiment  [11],  it 
has  been  subject  to  some  criticism.  One  of  the  weak  points  of  this  theory  is 
that  the  spanwise  wavenumber  of  the  three  dimensional  wave  was  chosen  to  fit 
the  experimental  data  rather  than  being  predicted  (see  Stuart  [12]). 

The  purpose  of  this  paper  is  to  follow  a  somewhat  different  approach  to  the 
cause  for  streamwise  vortices  in  a  turbulent  boundary  layer.  The  present 
hypothesis  sterns  from  what  is  known  as  the  Benney-Gustavsson  resonance  [13]  in 
r/hich  a  three-dimensional  disturbance  with  certain  wavenumbers  can  grow  to  a 
relatively  large  amplitude.  Benney  and  Gustavsson  showed  that  there  exist 
exact  resonances  in  various  linearized  laminar  stability  problems.  However, 
they  found  that  the  resonance  condition  is  only  approximately  satisfied  for 
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the  case  of  Blausius  flow.  The  effects  of  the  resonance  to  the  transition 
problem  is  thus  obscured  by  the  existence  of  an  unstable  Tollmien-Schlichting 
wave  and  the  fact  that  the  resonance  is  only  'close'. 

The  turbulent  pipe  flow  experiment  by  Morrison  and  Kronauer  [14]  indicates 
that  the  statistically  dominant  streamwise  fluctuations  have  a  v/avelike  cha¬ 
racter.  Further,  it  appears  that  a  linear  or  weakly  nonlinear  perturbation 
analysis  around  the  mean  velocity  might  be  applicable  to  the  turbulent  bound¬ 
ary  layer  problem.  Along  this  line,  Landahl  [15]  and,  subsequently.  Bark 
[16],  examined  the  implications  of  the  linear  theory  with  given  source  terms 
coming  from  an  assumed  model  for  the  Reynolds  stress.  A  particularly  interest¬ 
ing  result  is  that  Bark's  computed  energy  spectrum  for  the  streamwise  velocity 
fluctuations  displays  preferred  scales  in  the  spanwise  and  streamwise  direc¬ 
tions  and  in  time  which  are  in  good  agreement  with  Morrison  and  Kronauer 's 
measurements.  However,  the  bandwidth  of  the  computed  spectrum  is  a  few  orders 
of  magnitude  smaller  than  was  observed  experimentally. 

Bark  attributed  the  sharpness  of  his  computed  spectrum  to  the  crudeness  of  his 
model  for  the  Reynold's  stress.  However,  his  extremely  sharp  spectral  peak 
suggests  other  possibi 1 i ties  such  as  some  sort  of  resonance.  Bark's  method  of 
solving  the  linearized  equation  for  the  vertical  vorticity  implicitly  assumed 
that  there  is  no  resonance  such  as  suggested  by  Benney  and  Custavsson.  This 
raises  an  interesting  possibility  that  the  Benney-Gustavsson  resonance  becomes 
exact  if  one  replaces  the  Blasius  profile  with  turbulent  boundary  layer  pro¬ 
file.  By  examining  the  eigenvalues  of  the  Orr-Sonmerfeld  problem  and  the 
vertical  vorticity  equation,  this  was  found  to  be  the  case.  Details  will  be 
discussed  in  Section  2.  Since  this  resonance  becomes  exact  and  since  there  is 
no  unstable  mode  for  the  mean  velocity  profile,  the  effects  of  the  resonance 
could  be  more  pronounced  than  for  the  case  of  the  transition  problem.  Fur¬ 
thermore,  the  existence  of  this  resonance  seems  to  explain  the  observed  scales 
of  Morrison  and  Kronauer' s  experiment  and  the  extreme  sliarpness  of  Bark's  com¬ 
puted  spectral  peak. 
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i  th  this  explanation  for  the  preferred  scales  of  the  3-dimensional  disturb- 
ices,  we  then  examine  the  effects  of  the  nonlinear  interaction  on  this  reso- 
int  fundamental  mode  using  the  weakly  nonlinear  perturbation  technique.  It  is 
hown  in  Section  3  that,  as  an  integral  part  of  this  resonance,  there  is  a 
aan  secondary  flow  which  has  a  spanwise  wavelength  90  v/p /t^  and  whose 
elocities  exhibit  a  counter-rotating  streamwise  vortex  structure.  Further 
iscussions  are  made  in  Section  4. 
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n  Figure  7,  the  y-dependence  of  the  computed  vortex  structure  with  accompany- 
ng  low  speed  streaks  is  compared  with  the  vortex  structure  determined  with 
im’s  numerical  simulation  [17].  Only  the  shape  of  the  vortex  structure  is 
ompared  since  the  absolute  values  are  not  available  in  our  analysis  and  the 
elative  strengths  of  and  are  not  available  in  Kim’s  paper.  The 

imilarity  of  our  computed  vortex  shape  with  Kim's  simulation  supports  the 
elevance  of  the  resonance  mechanism  to  the  bursting  phenomena.  The  agreement 
etween  Kim's  simulation  results  and  Blackv/elder  and  Kaplan's  experimental 
esults  is  generally  good  if  the  slight  differences  are  attributed  to  differ- 
nces  between  temporal  and  spatial  averaging. 
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To  summarize,  the  induced  mean  flow  is  given  by 


/ 


'M  \  / 


/ 


U  (y,t)  cos(2B  z) 


Vm 


V,,(y,t)  cos(26^z) 

il  0 


V 


-1 


2BJ  3Vj^/3y  sin(28gZ) 


(3.22) 


where  U,^  is  related  to  ilf,|  by 

=  -  [1.^26^  .  (3.23) 

The  y-t  dependence  of  0^^  is  plotted  in  the  first  part  of  Figure  5.  As  for  the 
vertical  component  of  the  induced  mean  flow,  the  y-dependence  of  is  rather 
independent  of  time.  The  amplitude  of  Uj^  reaches  its  maximum  value  at 
t'*’  ■  150  and  then  decays  very  slowly  with  time.  In  order  to  interpret  the 
results  plotted  in  Figure  5,  it  is  helpful  to  plot  the  projection,  on  the 
(y,z)-plane,  of  the  streamlines.  Figure  6  shows  the  streamline  pattern  at 
t"*”  =  40.  Actually,  since  the  y-dependence  of  V.^  does  not  change  very  much 
with  time  either,  the  streamline  pattern  is  about  the  same  for  all  time.  The 
streamline  pattern  clearly  shows  the  counter-rotating  longitudinal  vortex 
structure  of  the  induced  mean  flow.  The  spanwise  wavelength,  of  the 
induced  mean  flow  is  90  which  compares  favorably  with  the  experimental  value 
of  100. 

At  z'*'  =  0,  the  counter-rotating  vortices  pump  the  low-speed  fluid  up  away  from 
the  wall  so  that  the  streamwise  component  of  the  induced  mean  flow  would  show 
a  momentum  defect.  These  low-speed  streaks  would  occur  every  90  wall  units  in 
the  spanwise  direction. 
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The  function  Vj^(t,y)  is  then  computed  by  taking  the  inverse  Laplace  transform 
where  the  integration  contour  in  the  complex  s-plane  should  lie  below  all  the 
poles  and  the  branch  cut  of  the  integrand.  Note  that  the  pole  associated  with 
the  source  term  is  located  on  the  branch  cut  since  R*  >  262/)(d  I  Further, 
the  singularities  in  the  expressions  for  the  coefficients  A^,  etc.,  are  remov¬ 
able.  Therefore,  the  contour  can  be  deformed  to  a  branch  cut  integral.  For 
the  integration  along  the  branch  cut  and  for  the  integration  over  y^,  we  have 
to  rely  on  numerical  evaluation.  Figure  5  (the  second  plot)  shows  the 
results.  Note  that  the  y-dependence  of  Vf;|(y,t)  does  not  change  very  much  with 
time.  The  amplitude  of  V|,^  reaches  its  maximum  value  at  t'*’  =  tu^/v  •<  40  and 
then  decays  very  slowly  compared  to  the  source  term's  decay  rate.  (Recall 
that  the  source  term  of  equation  (3.13)  decays  as  exp(-0.074  t'*’) ;  in  contrast, 
the  value  of  Vj,^  at  t"*"  =  200  is  still  about  1/3  of  its  peak.) 

To  obtain  the  horizontal  components  of  the  induced  mean  flow,  equation  (3.5) 

is  examined  up  to  the  order  of  Using  equation  (3.10)  for  and 

and  equation  (3.12)  for  the  equation  for  ^an  be  simplified  to 

U|t  ♦  "  I?)  -  ^  '‘1  "■  ''n  5in(26„z)  (3.20) 

As  for  the  computation  for  V[,|,  a  solution  is  sought  of  the  form 

n(i/2)  =  tlj^(y,t)  sin(2BoZ)  .  (3.21) 

The  resulting  equation  for  N,j,(y,t)  is  solved  again  as  an  initial  value  problem 
assuming  that  vanishes  initially.  Using  the  Laplace  transform  and  the 
Green's  function  technique,  the  computation  of  N.^  can  also  be  reduced  to  an 
evaluation  of  an  integral  over  y^  and  over  a  contour  in  the  frequency  domain. 
Once  N.^  is  computed,  the  horizontal  components  of  the  induced  mean  flow  are 
obtained  using  the  continuity  equation  and  the  definition  of  the  vertical 
vorti ci ty. 


«  ■■ 
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To  obtain  a  solution  Vf.^(y,t)  of  equation  (3.13),  which  vanishes  at  y  =  »,  we 
solve  equation  (3.13)  as  an  initial  value  problem  without  assuming  the  expo¬ 
nential  time  dependence.  Since  it  is  not  clear  which  initial  condition  for 
Vf.^(y,t)  is  an  appropriate  one,  we  shall  assume  that  it  vanishes  initially  and 
examine  the  y  and  t  dependences  of  the  resulting  solution. 


In  order  to  solve  the  initial  value  problem,  the  Laplace  transform  of  equation 
(3.13)  is  taken  to  give 


[is  (iL  -4sn  -  L  (li-  -4S2)^]  V  »  -4  BB*  * 

3y2  “  ^  3y2  “  ” 


3y  2(1)  j-i  s 


(3.17) 


where  V  is  the  Laplace  transform  of  V.^  defined  as 


^  -i  s t 

V,  (y.s)  =  /  dt  e  • 

0 


(3.18) 


The  above  equation  for  V.  is  solved  using  the  Green's  function  technique.  The 

M 

Green's  function  takes  the  following  form, 


-(4B2+iRs)^''^(y-y  )  -26(y-y^) 

A^e  °  +  B^e  ° 


for  y  >  y^ 


(3.19) 

-(4B2+iRs)^'^^(y-y  )  -2B(y-y.)  (4B2+iRs)^'^^(y-y  )  28(y-y„) 

A  e  °  +  B  e  °  +  C  e  °  +  D  e  ”  ° 


for  y  <  y 


where  y^  denotes  the  support  of  the  delta  function,  and  the  coefficients  A^, 
B^,  C_  and  D_  are  determined  from  the  boundary  condition  at  the  wall  and  the 
continuity  of  G,  G',  G"  and  the  jump  condition  on  G'"  at  y  =  y^.  Since  V.^ 
does  not  have  a  homogeneous  solution,  the  Green's  function  is  uniquely  defined 
by  these  conditions.  The  branch  cut  associated  with  the  square  root  is  chosen 
as  in  Figure  4. 
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where  uj  is  the  imaginary  part  of  hIq.  The  function  VQ2(y)  is  determined  from 


-  4»5)  -  2,a„S“  -  ^ 

ay2  *  3y^ 


=  -  2B2 


°o°l 


(3.14) 


As  long  as  (2(ijg)  is  not  an  eigenfrequency  of  the  0-S  problem  for  the  case  of 
a  =  Zoq  and  6=0,  the  above  equation  for  Vq2  has  a  solution. 


However,  the  equation  for  the  induced  mean  flow  requires  special  treatment. 
Although  the  source  term  in  the  equation  for  V.^  decays  with  time  as  exp(2u>jt), 
the  solution  Vj,|  cannot  have  the  same  time  dependence.  To  illustrate  this 
fact,  we  first  assume  that  Vf,]  has  the  same  exponential  time-dependence.  Then, 
equation  (3.13)  reduces  to 


12„, 


4  BB 


a2  62 


fa2+B2)2 
'•  0  o' 


3(V*o) 

8y 


(3.15) 


where  V,^  is  defined  by  V[^  =  Vj^(y)exp(2wjt).  For  simplicity  of  argument,  we 
assume  that  the  right-hand  side  of  the  above  equation  vanishes  for  large 
values  of  y.  The  asymptotic  solution  for  will  then  be  of  the  form. 


-2B  y  .(462+2<.  y 

V„  =  a  e  °  +  b  e 


(3.16) 


If  (2(OjR^+  4b2)  were  positive,  then  the  sign  of  the  square  root  could  be 
chosen  to  be  positive  so  that  would  vanish  at  infinity.  However,  the  argu¬ 
ment  of  the  square  root  at  the  resonant  point  is  negative. 
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where  and  8q  are  the  wavenumbers  at  the  direct  resonance.  Using  the  con¬ 
tinuity  equation  and  the  definition  of  N,  U^®^  and  can  be  obtained  from 
m(0)  as 


The  equation  for  is  obtained  from  equation  (3.4); 


v(0) 

3x2  ay 


+  2 


33 

3x3z3y 


(u(0),,(0)) 


33 


3z2  3y 


(.,,(0)|,(0)) 


(3.11) 


The  time  averages  of  and  are  dropped  since  U^*^^  and 

decay  exponentially  with  time. 


An  examination  of  the  source  terms  shows  that  the  particular  solution  for 
is  of  the  form 


/n\  2ia  x-2i(ij  t 

V  ^  =  Vj^(y,t)  coslZe^z)  +  (VQ2{y)  e  +  * 


(3.12) 


Vie  will  refer  to  the  first  term  as  the  induced  mean  flow  and  the  second  term 
as  the  induced  second  harmonic.  The  function  Vj^(y,t)  is  determined  as  a  solu¬ 
tion  of  the  following  equation 


1 


-  '•8?)')  * 
3y2  ° 


M 


4  BB* 


0  0 


t 


(3.13) 
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[' 


Ix^  IT 


^  A  N  +  el/2  u‘  ((UU-UU)  -  (UU-WW)) 


3z2  3x2  J 


+  e 


sk 


32 

3x3y 


(VW-VW) 


(3.5) 


3U 


+ 


el/2  il  = 
^  3y 


(3.6) 


N 


.  3U  3U 
■  ^  ■  3x 


(3.7) 


Equation  (3.5)  is  linear  when  e  =  0,  and  a  systematic  perturbation  scheme  is 
readily  developed  (i.e.,  N  =  +  £^^2  ^(1/2)  +  etc.).  The  zeroth 

order  term  in  such  an  expansion  gives 

*  .k  .  .  * 

fnl  «  *  -  “I® 

tj  =  B  Tig(y)  sin  e^z  e  +  B  r\^(y}  sin  B^z  e  °  °  (3.8) 


where  *  denotes  the  complex  conjugate.  The  function  n^  is  determined  as  a 
solution  to  the  eigenvalue  problem 


[ia^(i]-c) 


'ii 

■3y2 


=  0 


n^(0)  = 


)  =  0 


(3.9) 
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3.  MEAN  FLOW  INDUCED  BY  THE  RESONATING  FUNDAMENTAL  MODES 

From  equations  (2.15)  and  (2.16),  it  is  apparent  that  if  a,  e  and  m  correspond 
to  a  direct  resonance,  then  so  do  a,  -6  and  «>.  So,  there  are  at  least  two 
resonant  modes.  The  nonlinear  theoretical  consequences  of  two  or  more  reso¬ 
nant  modes  were  analyzed  by  Benney  and  Gustavsson.  We  shall  rely  on  their 
formalism  to  compute  the  mean  flow  induced  by  these  resonant  modes. 


They  first  observed  that  equations  (2.6)  and  (2.7)  imply  that,  for  the  case  of 
two  resonant  modes,  the  dominant  terms  in  the  higher-order  expansion  for  v  and 
n  are  of  the  form 


V  ~  (l  +  e^t**  +  ...)  e 


iax+i6z-i(ut 


(3.1) 


( t  +  t^  +  . . . )  e 


iax+i6z-i(i)t 


(3.2) 


where  e  is  the  nonlinear  parameter  related  to  the  small  amplitude  of  fluctua¬ 
tion.  Thus,  they  argue,  the  appropriate  time  scale  is  and  the  perturba¬ 
tion  must  be  rescaled  in  the  form 


V  =  V,  (u,w,n)  =  (u,H,N) 


(3.3) 


The  nonlinear  perturbation  equations  then  become 

=  ‘  It  -  ' 


[‘ 


3X^ 


3x 


33 


3x23y 


(UU-UU]  -  2 


3^  /'...I  i.i.-i  3^ 


3x3z3y 


(UW-UW)  - 


Sz^ay 


(WU-WU) 


+  e 


1/2 


(IL  +  -  3^)  (9  (UV-UV)  +  3  (WV-WV)) 

L  3x2  3z2  3v2  ^ 


+  e  (ii-  +  — )  (VV-VV)  =  0 

L  3x2  3z2  3y  J 


(3.4) 


-11- 


Dynamics  Technology 
DT-8154-06 


density  as  a  function  of  and  either  a  or  B.  Figure  3  shows  their  contour 
plots  of  3?(a^,(ij^)  and  where  these  quantities  are  related  to  the 

power  spectrum  as 

=  a*' tii*'  (2.27) 


^(B'^.o)'^)  =  P^(B'^.uj'^)  (2.28) 

These  spectral  functions  plotted  in  logarithm  scales  permit  one  to  make  a 
visual  comparison  of  the  relative  power  in  various  wavenumber/frequency  bands 
(the  power  is  the  contour  level  times  the  area  as  seen  in  the  figure). 

The  particular  spectra  measured  at  y'*'  =  14.8  were  chosen  because  this  location 
is  the  closest  of  their  data  to  the  peak  of  the  resonant  free  mode  of  the 
vertical  vorticity  equation  (see  Figure  1).  The  predicted  position  of  the 
spectral  peak  based  on  the  direct  resonance  is  denoted  by  the  intersection  of 
the  two  dashed  lines  in  each  spectrum.  The  prediction  is  in  good  agreement 
with  the  experimental  data. 


A  ^  A 

e  /  TipU'  V  dy 
I  Tip  rip  dy 
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(2.23) 


We  examined  whether  such  a  resonance  occurs  when  u  is  taken  to  be  a  turbulent 
boundary  layer's  mean  profile.  For  this  purpose,  we  have  to  rely  on  numerical 
solutions  of  the  eigenvalue  problems  associated  with  equations  (2.15)  and 
(2.18).  If  two  eigenvalues  coincide,  then  resonant  growth  occurs.  One  such 
resonance  was  found  when  the  wavenumbers  (expressed  in  'wall  units'; 
a'*’  =  av/u_,  h  Bv/u  ,  where  u  e  /t  /p  )are 

‘  T  T  W  ' 


0.0093 


B  =  0.035 


(2.24) 


(2.25) 


The  location  of  this  resonance  in  this  dimensionless  wavenumber  appears  to  be 
independent  of  the  Reynolds  number  (vre  examined  the  Reynolds  number  dependence 
for  1000  <  R*  <  15,000).  At  this  resonant  point  the  eigenfrequencies  of  the 
0-S  and  the  V-V  equations  are  the  same  and  given  by 


0,"^  =  ^  =  0.090  -  i0.037 


(2.26) 


This  resonance  was  found  numerically  by  locating  the  intersection  of  the  zero 
lines  of  the  real  and  the  imaginary  parts  of  [c(a,B;R*)  -  c'(a,B;R*)]  in  the 
a-B  plane.  The  boundary  layer  velocity  profile  used  for  these  computations 
incorporated  the  well-known  Law  of  the  Wall  and  Law  of  the  Wake  along  with  a 
sublayer  "patch"  based  on  van  Driest's  damping  factor.  This  sublayer  patch 
results  in  a  velocity  profile  that  is  in  good  agreement  with  experiment  over 
the  entire  boundary  layer.  The  normal  coordinate  dependence  of  the  0-S  eigen- 
mode  and  V-V  eigenmode  are  plotted  in  Figures  1  and  2,  respectively. 

Because  of  the  potentially  large  amplitude  associated  with  the  secular 
behavior,  the  above  values  of  a'^,  B'*'  and  (real  part  of  oj'*')  should  corres¬ 
pond  to  the  position  of  a  local  peak  in  the  power  spectral  density  of  the 
horizontal  velocity.  Experimental  evidence  for  the  correctness  of  this  pre¬ 
diction  can  be  found  in  Morrison  and  Kronaucr's  experiment  [14].  Measuring 
the  fluctuation  of  the  streamwise  velocity,  they  obtained  the  power  spectral 


J 
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subject  to  the  boundary  conditions;  v=0,^=0,n=0aty  =  0,  «.  Here, 
c  denotes  the  phase  velocity  defined  by  c  =  w/ct.  Equation  (2.15)  is  the  Orr- 
Sommerfeld  (0-S)  equation,  and  (2.16)  will  be  referred  to  as  the  vertical  vor- 
ticity  (V-V)  equation.  The  0-S  eigenvalue  problem  leads  to  a  set  of 
eigenvalues 


c  =  c(a,  B;  R*]  . 


(2.17) 


The  V-V  equation  is  usually  solved  as  a  forced  response.  This  approach 
implicitly  assumes  that  the  eigenfrequencies  of  the  0-S  equation  do  not  match 
any  of  the  free  mode  eigenfrequencies  for  the  V-V  equation,  that  is,  the 
eigenfrequencies  of  the  problem 


[io(u-c)  -  ^  (—  -  a2-B21]  =  0  ;  nr  =  0  at  y  =  0,  «  . 

3y2  ^ 


Denoting  the  set  of  V-V  free-mode  eigenfrequencies  by 


(2.18) 


c  =  c  (o»8;R*)  > 


the  condition  for  direct  resonance  can  be  written  as 


c(a,B;R*)  =  c'(o,B;R*)  . 


(2.19) 


(2.20) 


If  such  a  resonant  condition  is  satisfied,  the  solution  to  the  V-V  equation 
behaves  as 


n  ~  at  np(y)  e' 


■iu)t+iax+iBz 


while 


V  ~  v(y)  e 


-iut+iax+iBz 


(2.21) 


(2.22) 


Mote  that  the  secular  solution  for  n  initially  grows  linearly  with  time  so 
that  relatively  large  horizontal  motions  might  be  induced.  Here,  the  parame¬ 
ter  'a'  which  determines  the  initial  growth  rate  of  the  secular  term  is 
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Si  =  (uu  -  uu)  +  ly  (vu  -  vu)  +  (wu  -  wu)  (2.8) 

^2  =  I7  (iJV  -  uv)  +  ly  (vv  -  vv)  +  |y  [m  -  wv)  (2.9) 

S3  =  -jy  (uw  -  uw]  +  _  (vw  -  vw)  +  _  (vM  -  ww)  .  (2.10) 


The  horizontal  components  of  the  velocity  perturbation  are  related  to  v  and  n 
by  the  equations 


au  .  av;  _  av 

■5x  "Sz  ■  '  ly 


(2.11) 


au 

Jz 


(2.12) 


Ue  next  linearize  equation  (2.3)  and  seek  normal  mode  solutions  of  the  form 


V  =  v(y) 


(2.13) 


n  =  n(y) 


giax+iBz-iut 


(2.14) 


where  o  and  B  are  wavenumber  components  in  the  streamwise  and  spanwise  direc¬ 
tion  and  to  is  the  wave's  frequency.  The  following  equations  are  then  obtained 


[ia(u-c) 


a2-B2)  -  iaH"  -  ^  (^  -  a2-B2)]  ‘v  =  0 

ay2 


(2.15) 


1 


J 


[ia(u-c)  -  ^  -  o2-82)]  n  =  -  iBu'  v 

*  ay2 


(2.16) 
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resonance  by  which  disturbances  with  certain  scales  can  grow,  at  least  momen¬ 
tarily,  by  subtracting  energy  from  the  mean  flows.  Since  the  eigenmodes  whose 
peak  values  are  near  the  turbulent  boundary  layer's  sublayer  decay  rather 
rapidly,  such  a  resonance  will  only  be  effective  if  it  occure  in  the  lowest 
possible  order  of  the  expansion  in  order  to  reverse  the  energy  loss  even 
momentarily. 

The  direct  resonance  suggested  by  Benney  and  Gustavsson  does  arise  ir  the 
linearized  equations  of  (2.3)  and  (2.4).  Another  interesting  property  of  the 
direct  resonance  is  that  at  a  given  Reynolds  number,  it  occurs  only  at  dis¬ 
crete  locations  in  wavenumber  space  (i.e.,  only  those  components  with  certain 
combinations  of  wavelength  and  obliquity  are  resonant).  Therefore,  this  reso¬ 
nance  seems  to  be  an  ideal  candidate  for  a  mechanism  which  allows  disturbances 
with  certain  scales  to  be  energetically  dominant. 

To  examine  the  possibility  of  direct  resonance  in  a  turbulent  boundary  layer, 
the  small  disturbance  equations  are  made  dimensionless  by  utilizing  the  bound¬ 
ary  layer's  displacement  thickness,  6*,  and  the  freestream  velocity  u^.  In 
the  non-dimensional  form  of  the  equations,  we  set  u^i5*/v  s  R*.  As  is 
customary  for  the  boundary  layer  problem,  x,  y  and  z  are  taken  to  denote  the 
streamwise,  the  normal  and  the  spanwise  coordinates,  respectively.  The  veloc¬ 
ity  components  will  be  denoted  by  u  ,  v  and  w.  Elimination  of  the  pressure 
from  (2.3)  gives 


(It  ^  Iv)  Av 


■57  J 


7."  , 

Ix 


9^  S3 
9z9y^ 


a2 


(2.6) 


u 


9 

■57 


]  n  +  u' 


9v 

Jz 


+  e 


9  S 1 


9S3 

TT 


1  =1 


An 


(2.7) 


flere  the  prime  denotes  differentiation  with  respect  to  y,  and  n  is  the  ver¬ 
tical  vorticity  ('i  "  §7)  •  e  is  a  nonlinear  parameter  related  to  the 

small  amplitude  of  the  fluctuations,  and  Sj  are  defined  by 
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2.  DIRECT  RESONANCE  IN  A  TURBULENT  DOUMDARY  LAYER 

The  velocity  and  pressure  field  in  a  turbulent  boundary  layer  are  split  into  a 
mean  and  a  fluctuating  part  according  to 


U.  =  u.  +  u. 
1  1  ^ 


P  =  p  +  p 


(2.1) 


(2.2) 


where  the  overbar  denotes  the  time  average.  Then,  the  equations  governing  the 
fluctuating  part  can  be  obtained  by  substituting  (2.1)  and  (2.2)  into  the 
Navier-Stokes  equations  and  subtracting  the  mean  part; 


^  .  U.7.U.  .  U.V.U.  = 


-v.p  +  vau^.  +  V.  (Uju^-u.u,.) 


(2.3) 


V.u^  =  0  . 


(2.4) 


Here,  p  is  the  density  of  the  fluid  and  v  is  the  kinematic  viscosity,  both 
taken  as  constant. 


The  mean  flow  will  be  assumed  to  be  parallel  such  that 


=  (u(y),  0,  0)  . 


(2.5) 


The  mean  velocity  distribution  u(y)  for  the  case  of  a  turbulent  boundary  layer 
will  be  considered  to  be  known.  The  key  assumption  in  our  analysis  of  (2.3) 
is  that  the  effects  of  nonlinear  terms  are  important  only  intermittently.  The 
streamwise  vortices  and  the  accompanying  low-speed  streaks  which  seem  to  pro¬ 
vide  the  set-up  for  a  burst  are  assumed  to  be  weakly  nonlinear  phenomena. 


With  these  assumptions,  equations  (2.3)  and  (2.4)  are  analyzed  using  perturba¬ 
tion  techniques.  It  is  v/ell  known  that  all  the  eigenvalues  of  the  linearized 
equations  of  (2.3)  and  (2.4)  result  in  decaying  solutions,  i.e.,  all  disturb¬ 
ances  lose  energy  to  the  mean  flow.  However,  there  might  be  a  possible 
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4.  CONCLUDING  REMARKS 

It  is  natural  to  be  concerned  with  the  robustness  of  the  direct  resonance,  and 
the  effects  on  the  resonant  condition  of  different  boundary  conditions,  mean 
profiles,  and  with  eddy  viscosity  terms  in  the  equation  were  examined. 
Including  eddy  viscosity  terms  does  not  affect  the  existence  of  a  direct 
resonance  and  the  computed  induced  mean  flow  exhibits  features  similar  to 
those  described  in  Section  3.  With  the  mean  profile  measured  by  Reischman  and 
Tiderman  [18]  for  a  polymer  injected  turbulent  boundary  layer,  the  direct 
resonance  was  still  found  to  remain.  The  spanwise  wavenumber  at  the  resonance 
is,  however,  smaller  than  that  for  the  'universal'  mean  profile,  which  agrees 
with  experimental  evidence  of  the  wider  streak  spacing  for  a  polymer  added 
flow.  We  also  showed  that  direct  resonance  should  be  virtually  unaffected 
even  if  the  rigid  wall  were  replaced  with  a  compliant  wall. 

As  will  be  discussed,  there  are,  of  course,  several  aspects  of  the 
theory  that  remain  incomplete.  However,  this  study  provides  strong  evidence 
for  the  relevance  of  direct  resonance  to  the  streamwise  vortices  observed  in  a 
turbulent  boundary  layer.  In  terms  of  the  energy  exchange  between  the  mean 
flow  and  the  fluctuations,  direct  resonance  can  be  interpreted  as  follows: 
the  source  term  in  the  linearized  vertical  vorticity  equation  (2.16)  reflects 
the  production  of  horizontal  disturbance  energy  by  the  action  of  the  mean 
field,  u,  against  the  uv-component  of  the  Reynold's  stresses.  When  the  wave- 
numbers  and  the  frequency  associated  with  an  0-S  eigenmode  correspond  to  those 
of  a  free  mode  of  the  vertical  vorticity,  the  production  of  energy  through  the 
uv-Reynolds  stress  becomes  more  efficient.  This  situation  is  similar  to  the 
harmonic  oscillator  driven  by  broadband  noise  about  its  natural  frequency. 

There  are  two  main  pieces  of  evidence  for  the  relevance  of  such  a  mechanism  to 
the  bursting  process.  One  is  that  the  wavenumbers  and  frequency  at  resonance 
are  close  to  the  values  associated  with  the  most  intense  v/aves  measured  by 
Morrison  and  Kronauer  near  the  sublayer  boundary.  The  other  is  that  the 
secondary  mean  flow  induced  by  this  resonant  fundamental  mode  contains  a 
streamwise  vortex  structure.  The  theoretical  shape  of  the  vortices  and  the 
spacing  of  the  accompanying  low-speed  streaks  are  comparable  to  experimental 
findings. 
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One  of  the  problems  which  has  not  been  addressed  is  a  comparison  of  the  reso¬ 
nant  point  and  other  points  in  wavenumber  space.  In  order  to  extend  the 
theory  to  predict  the  power  spectral  distribution,  this  issue  has  to  be 
addressed.  A  difficulty  associated  with  this  extension  is  that  although  the 
resonant  mode  initially  grows  with  time,  the  exponentially  decaying  factor 
eventually  dominates  and  the  mode  decays.  Therefore,  without  a  fresh  supply 
of  disturbances,  this  resonance  cannot  sustain  itself.  In  this  respect,  a 
potentially  important  property  of  this  resonant  mode  seems  to  be  the  inflec¬ 
tional  instability.  As  a  test,  we  considered  a  one  parameter  family  of  mean 
profiles  which  starts  with  the  universal  mean  profile  at  s  =  0  and  ends  with 
the  inflectional  profile  found  by  Blackwelder  and  Kaplan  [6]  at  s  =  1,  [i.e., 
Ugly)  =  (1-s)  u(y)  +  s  Ug|^,  where  Ugj/  denotes  the  Blackwelder  and  Kaplan  pro¬ 
file).  As  s  increases,  the  damping  factor  associated  with  the  Orr-Sommerfel d 
mode  becomes  smaller  and  near  s  =  0.7  that  mode  becomes  unstable.  As  it 
becomes  unstable,  it  can  no  longer  resonate  with  the  vertical  vorticity  mode 
since  it  can  be  shown  analytically  that  all  the  free  vertical  vorticity  modes 
are  damped.  Although  at  s  =  0  the  resonant  Orr-Sommerfeld  mode  is  not  the 
least  damped  mode  (actually  it  is  the  second  least  damped  mode;  the  eigen- 
frequency  of  the  least  damped  mode  being  u'*'  =  0.13  +  i0.035  at  o'*'  and  given 
by  equations  (2.24)  and  (2.25)),  as  s  increases  this  order  in  damping  factor 
changes  and  only  the  resonant  mode  becomes  unstable  at  s  =  1.  These  two  pro¬ 
perties  of  the  second  mode  (when  ordered  by  the  imaginary  part  of  the  eigen- 
frequency),  i.e.,  the  direct  resonance  and  the  inflectional  instability,  might 
be  the  reasons  why  it  plays  an  important  role  even  if  it  has  a  slightly  larger 
damping  factor.  However,  to  estimate  these  effects  quantitatively  appears  to 
require  a  more  systematic  account  of  the  several  lowest  modes  with  varying 
wavenumbers. 

Another  problem  which  has  not  been  considered  concerns  the  origin  and  math¬ 
ematical  description  of  the  sweep.  This  might  involve  a  very  complicated 
interaction  between  the  outer  flow  and  the  near  v/all  flow.  In  this  regard,  it 
seems  appropriate  to  mention  that  the  induced  mean  flow  in  Section  3  was 
represented  by  an  integral  along  a  branch  cut  in  the  complex  frequency  plane. 
This  branch  cut  is  what  Mack  [19]  interpreted  as  the  continuous  spectrum.  For 
a  frequency  in  this  continuous  spectrum,  the  corresponding  eigensolution  does 
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not  decay  exponentially  but  becomes  oscillatory  as  y  approaches  infinity.  This 
raises  the  possibility  that  the  interaction  between  the  inner  and  the  outer 
flows  could  be  described  through  the  induced  mean  flow.  Recall  that  the 
initial  condition  for  the  induced  mean  flow  was  chosen  to  be  zero  in  Section  3 
since  we  do  not  know  what  an  appropriate  condition  should  be.  Perhaps  by 
choosing  an  initial  condition  appropriate  to  a  finite  wavepacket  rather  than 
an  infinite  plane  wave,  the  time  development  of  the  induced  mean  flow  would 
exhibit  characteristics  of  the  sweep  in  addition  to  the  vortex  structure. 
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Figure  4.  Branch  cut  of  (46^  +  (4^2  +  iRS)i'2  is  defined  by 

R^'2  r  exp  i  (0/2  +  11/4)  where  0  varies  from  (-3it/2)  to  (ir/2). 
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